Optimal investment under multiple defaults risk: 
a BSDE-decomposition approach 

Ying JIAO* Idris KHARROUBit Huyen PHAM* 
January 13, 2013 

Abstract 

We study an optimal investment problem under contagion risk in a financial model 
subject to multiple jumps and defaults. The global market information is formulated 
as progressive enlargement of a default-free Brownian filtration, and the dependence 
of default times is modelled by a conditional density hypothesis. In this Ito-jump 
process model, we give a decomposition of the corresponding stochastic control problem 
into stochastic control problems in the default-free filtration, which are determined 
in a backward induction. The dynamic programming method leads to a backward 
recursive system of quadratic Backward Stochastic Differential Equations (BSDEs) in 
Brownian filtration, and our main result is to prove under fairly general conditions the 
existence and uniqueness of a solution to this system, which characterizes explicitly the 
value function and optimal strategics to the optimal investment problem. We illustrate 
our solutions approach with some numerical tests emphasizing the impact of default 
intensities, loss or gain at defaults, and correlation between assets. Beyond the financial 
problem, our decomposition approach provides a new perspective for solving quadratic 
BSDEs with finite number of jumps. 
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1 Introduction 



In this paper, we address an investment problem in an assets portfolio subject to defaults 
and contagion risk, which is a major issue for risk management in financial crisis period. We 
consider multiple default events corresponding for example to the defaults of multi credit 
names or to counterparty defaults, and contagion effects meaning that defaults on some 
assets may induce loss or gain on the other assets. One usually formulates the default-free 
assets price process as an Ito process governed by some Brownian motion W, and jumps are 
introduced at random default times, associated to a marked point process /i. The optimal 
investment problem in this incomplete market framework may be then studied by stochastic 
control and dynamic programming methods in the global filtration G generated by W and 

This leads in principle to Hamilton-Jacobi-Bellman integrodifferential equations in a 
Markovian framework, and more generally to Backward Stochastic Differential Equations 
(BSDEs) with jumps, and the derivation relies on a martingale representation under G 
with respect to W and fx, which holds under intensity hypothesis on the defaults, and 
the so-called immersion property (or (H) hypothesis). Such an approach was used in the 
recent papers [1], [13] in the single default case, and in [7] for the multiple defaults case. 
For exponential utility criterion, the solution to the optimal investment problem is then 
characterized through a quadratic BSDE with jumps, whose existence is proved under a 
boundedness condition on the portfolio constraint set. 

We revisit and extend the optimal investment problem in this multiple defaults context 
by using an approach initiated in [8] in the single default time case, and further developed 
in [H] in the multiple defaults with random marks case. By viewing the global filtration G 
as a progressive enlargement of filtrations of the default-free filtration F generated by the 
Brownian motion W, with the default filtration generated by the random times and jumps, 
the basic idea is to split the global optimal investment problem, into sub-control problems 
in the reference filtration F and corresponding to optimal investment problems in default- 
free markets between two default times. More precisely, we derive a backward recursive 
decomposition by starting from the optimal investment problem when all defaults occurred, 
and then going back to the initial optimal investment problem before any default. The main 
point is to connect this family of stochastic control problems in the F-filtration, and this 
is achieved by assuming the existence of a conditional density on the default times given 
the default-free information F. Such a density hypothesis, which is standard in the theory 
of enlargement of filtrations, was recently introduced in [3], [5] for credit risk analysis, and 
may be seen as an extension of the usual intensity hypothesis. 

This F-decomposition approach allows us furthermore to formulate an optimal invest- 
ment problem where the portfolio constraint set can be updated after each default time, 
depending possibly on the past defaults, which is financially relevant. This extends the 
global approach formulation where the portfolio set has to be fixed at the beginning. Next, 
for exponential utility function criterion, we apply dynamic programming method to each 
optimal investment problems in the F-filtration. We then get rid of the jump terms arising 
in the dynamic programming in the G-filtration, and are led instead to a backward recursive 
system of quadratic BSDEs in Brownian filtration with a nonstandard exponential term. 
Our main result is to prove under fairly general conditions (without assuming in particular 
a boundedness condition on the portfolio constraint set) the existence and uniqueness of a 
solution to this system of BSDEs. Existence is showed by induction, based on Kobylanski 
results [12] together with approximating sequences for dealing with the exponential term 
and unbounded portfolio, suitable uniform estimates and comparison results for getting the 
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convergence. Uniqueness is obtained by verification arguments for relating the solution of 
these BSDEs to the value functions of the F-control problems, and uses BMO-martingale 
tools. Moreover, an interesting feature of our decomposition is to provide a nice charac- 
terization of the optimal trading strategy between two default times, and to emphasize the 
impact of defaults and jumps in the portfolio investment. We also illustrate numerically 
these results in a simple two defaultable assets model, where each asset is subject to its own 
default and also to its counterpart. Finally, we mention that beyond the optimal investment 
problem, the F-decomposition approach provides a new perspective for solving (quadratic) 
BSDEs with finite number of jumps, see the recent paper 

The outline of this paper is organized as follows. In Section 2, we present the multiple 
defaults model where the assets price process is written as a change of regimes model with 
jumps related to the default times and random marks. Section 3 formulates the optimal 
investment problem, and gives the decomposition of the corresponding stochastic control 
problem. Section 4 is devoted to the derivation by dynamic programming method of the 
sub-control problems in terms of a recursive system of BSDEs, and to the existence and 
characterization results of this system for the optimal investment problem. Finally, we 
provide in Section 5 some numerical experiments for illustrating our solutions approach in 
a simple two-defaultable assets model. 

2 Multiple defaults model 
2.1 Market information setup 

We fix a probability space {U,Q,F), equipped with a reference filtration F = {J-'t)t>o sat- 
isfying the usual conditions, and representing the default-free information on the market. 
Let T = (ri, • • • ,r„) be a vector of n random times, representing multiple default times, 
and L = (Li, ■ ■ ■ ,Ln) be a vector of n marks associated to default times, Lj being an G- 
measurable random variable taking values in some Polish space £" C M^', and representing 
for example the loss given default at time Tj. The global market information is given by 
the default-free information together with the observation of the default times and their 
associated marks when they occur. It is then formalized by the progressive enlargement of 
filtration G = F V . . . V B", where D'' = {V^)t>o, = i>t+ , i>t = cr{Lklrk<s, s < t) y 
(t(1t-j.<s, s < t), k = 1, . . . ,n. In other words, G = {Qt)t>o is the smallest right-continuous 
filtration containing F such that for any /c = 1, . . . , n, is a G-stopping time and is 
^T-fc-measurable. 

For simplicity of presentation, we shall assume in the rest of this paper that the default 
times are ordered, i.e. ti < . . . < r„, and so valued in A„ on {t„ < oo} where 



On one hand, this means that we do not distinguish specific credit names, and only observe 
the successive default times, which is relevant in practice for classical portfolio derivatives 
like basket default swaps. On the other hand, we may notice that the general non-ordered 
multiple random times case for (ti, . . . , r^) (together with marks (Li, . . . , can be de- 
rived from the successive random times case by considering suitable auxiliary marks. In- 
deed, denote by fi < . . . < the corresponding ordered times, and by Lk the index mark 
valued in {1, . . . , n} so that = t^^ for k = 1, . . . ,n. Then, it is clear that the progressive 
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enlargement of filtration of F with the successive random times (fi, . . . ,-?■„) together with 
the marks (ii, , . . . , i„, L^^) leads to the filtration G. 

We introduce some notations used throughout the paper. For any (^i, . . . ,9n) G A„, 
(4, . . . ,4) e E^, we denote by 6> = (^i,. . .,9„), £ = {h, . . . and Ok = {Oi, . . .,6k), 4 
= {£i, . . . , for /c = 0, . . . , n, with the convention Oq = = 0. We also denote by = 
(n, . . . , Tfc), and Lfc = (Li, . . . , L^). For t > 0, the set ilj denotes the event 

■■= {rk<t<Tk+i}, 

(with = {t < Ti}, il" = {r„ < t}) and represents the scenario where k defaults occur 
before time t. We call as the A;-default scenario at time t. We define similarly = 
{^k <t< Tk+i}. Notice that for fixed t, the family {^t)k=o,...,n (resp. {^^-)k=o,...,n) forms 
a partition of 0. We denote by ■p(F) the cr- algebra of F-predictable measurable subsets on 
R+ X Q, and by Vvi^^, E^) the set of indexed F-predictable processes Z^{., .), i.e. s.t. the 
map {t,Lo,6k,£k) Z^{L0,9k,£k) is 'P(F) (g) 5(Afc) (g) ;6(£j'^)-measurable. We also denote by 
OwiA'', E^) the set of indexed F-adapted processes Z .), i.e. such that for all t > 0, the 
map {u},6k,ik) — > Z^{uj,Ok,ik) is Tt ^ B{Ak) (g) ;S(E^)-measurable. 

We recall from [14:\ Lem2.1] or [9l Lem4.1] the key decomposition of any G-adapted 
(resp. G-predictable) process Z = (^t)t>o in the form 

n n 

Zt = Y,l^,Z^{Tk,Lk), (resp. = J^l^._Zf(Tfc,i:fc)), i > 0, 

fc=0 fc=0 * 

where Z^ lies in Ow{Ak,E^) (resp. V^iAk^E^)). 

As in [5] and [H], we now suppose the existence of a conditional joint density for (r, L) 
with respect to the filtration F. 

Density Hypothesis. There exists a G 0^{lS.n,E^) such that for any bounded Borel 
function / on A" x E^, and t > 0: 

nf{r,L)\Tt] = [ f{e,l)at{e,l)der,{dt) a.s., (2.1) 
where d9 = dOi . . . dO^ is the Lebesgue measure on M" , and 'q{d€) is a Borel measure on E"^ 

n-l 

in the form ri[dt) = r]i{dii) ?7fc+i(£fc, di^+i), with rji a nonnegative Borel measure on E 

k=l 

and r}k+i{tk,dlk+i) a nonnegative transition kernel on E^ x E. 

Remark 2.1 From the condition (12. ip . we see that r admits a conditional (w.r.t. F) 
density with respect to the Lebesgue measure given by a'^{6) = J a{0, i)r]{d£). This implies 
in particular that the default times are totally inacessible with respect to the default-free 
information, which is consistent with the financial modelling that the default events should 
arrive by surprise, and cannot be read or predicted from the reference market observation. 
This joint density condition w.r.t. the Lebesgue measure also implies that the default times 
cannot occur simultaneously, i.e. Ti tj, i ^ j, a.s., which is a standard hypothesis in the 
modelling of multiple defaults. Moreover, by considering a conditional density, and thus a 
time-dependence of the martingale density process {at{6, i))t>o, we embed the relevant case 
in practice when the default times are not independent of the reference market information 
F. Compared to the classical default intensity processes for successive defaults in the top- 
down modelling approach, the conditional density provides more and necessary information 
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for analyzing the impact of default events. Further detailed discussion and some explicit 
models for density of ordered random times are given in [S]. 

On the other hand, the condition ()2.ip implies that the family of marks L admits a con- 
ditional (w.r.t. F) density with respect to the measure r]{d£) given by a^{£) = f a{6,£)d6. 
This general density hypothesis (j2.ip embeds several models of interest in applications. In 
the case where a is separable in the form a{0,£) = {6)a^ {£), this means that the ran- 
dom times and marks are independent given Tt ■ The particular case of nonrandom constant 
mark = ik is obtained by taking Dirac measure rjk = The case of i.i.d. marks L^, 
/c = 0, . . . , ?i is included by taking a^{Z) separable in and 77 as a product measure. We 
can also recover a density modelling of ordered default times (as in the top-down approach) 
from a density model of the nonordered defaults (as in the bottom-up approach). Indeed, 
let T = (ti, • • • , Tn) be a family of nonordered default times having a density a"^, and denote 
by T = (fi, . . . , f„), L = {li, . . . ,Ln) the associated ranked default times and index marks. 
By using statistics order, we then see that (t, l) satisfy the density hypothesis with 

a{6l, . . . ,6n,k, ■ ■ ■ ,in) = ^ CH.'^ {0a{l)i ■ ■ ■ ,Gcr{n))'^{{ii,...,i^) = {cr{l),...,cr{n))} 

for (^1, . . . , On) € An, i = {ii, . . . ,in) € E = {!,..., n}, where S.„ denotes the set of all 
permutations a = {a{l), . . . ,a{n)) of E, and with i]{M) = ^^^^^ '^f=(T; i]k+ii^k,d^) = 

2.2 Assets and credit derivatives model 

We consider a portfolio of d assets with value process defined by a d-dimensional G-adapted 
process S. This process has the following decomposed form 

n 

St = ^l^,S^{Tk,Lk), (2.2) 

k=0 

where S^{6k, ik)-, ^fc = (^1, ■ ■ ■ , ^fc) S A^, = {£1, . . . , Ik) S E^ , is an indexed process in 
OF(Afc, E^), valued in M^^, representing the assets value in the /c-default scenario, given the 
past default events = Ok and the marks at default = £k- Notice that St is equal 
to the value S^ only on the set Q^l, that is, only for < t < t^+i- We suppose that the 
dynamics of the indexed process S^ is given by 

d5f(0fc,4) = s1{ek,ik)*{bt{0k.e-k)dt + <j'i{ek,ik)dWt), tye^, (2.3) 

where is a m-dimensional (P, F)-Brownian motion, m > d, and o"'^ are indexed pro- 
cesses in Vwi^k, E''), valued respectively in and R'^^"^. Here, for 2; = {xi, . . . , Xd)' G M'^, 
and y = (yi, . . . , y^)' in M'^^'^, the expression x * y denotes the vector {xiyi, . . . , Xdyd)' in 
j^dxg rjij^g model ()2.2p - ()2.3p can be viewed as an assets model with change of regimes after 
each default event, with coefficients , depending on the past default times and marks. 
We make the usual no-arbitrage assumption that there exists an indexed risk premium 
process A'^ € Vw{Ak,E^) s.t. for all (0^,4) € x E'' . 

aHOk,ik)Xt{Ok,h) = bU0k,h), t>0. (2.4) 

Moreover, in this contagion risk model, each default time may induce a jump in the assets 
portfolio. This is formalized by considering a family of indexed processes 7'^, /c = 0, . . . , n — 
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1, in Vf{A'',E'',E), and valued in [-1,00)"^. For (0^,4) G A'^ x E^, and 4+i G E, 
JtiO^jikj^k+i) represents the relative vector jump size on the d assets at time t = O^+i > 
9k with a mark ik+i, given the past default events {rk^Lk) = (^fc^-^fe)- In other words, we 
have: 

S'^+UOk+i^lk+i) = (0fc,4)*(ld + 7e\+,(0fc,4,4+i)), (2.5) 

' fc+ 1 

where we denote \d as the vector in with all components equal to 1. 

Remark 2.2 In this defaults market model, some assets may not be traded anymore af- 
ter default times, which means that their relative jump size is equal to —1. For k = 
0, ...,n, {dk,ik) € Afc X E'^, denote by d^{6k,^k) the number of assets among the d- 
assets which cannot be traded anymore after k defaults, so that we can assume w.l.o.g. 
6'=(0fc,4) = (fo'(0fc,4) 0), a^{ek,tk) = (a'(0fc,4) 0), 7'(^fc,4,^) = (7^(0fc,4,^) 0), 
where 6'^(0fc,4)) ^^(^fc5 4)) l^i'^kj^ki^) are F-predictable processes valued respectively in 
^d^OkM^ I^J'=(0fcA)xm^ ]Krf*(0feA) with d^{ek,ik) = d-d''{ek,ik), the number of remaining 
tradable assets. Either dJ'{Ok,ik) = 0, and so a^{6k,tk) = 0, b''{Ok,ik) = 0, j'^iOkjikJ) 
= 0, in which case (12. 4p is trivially satisfied, or d^{Ok-,ik) ^ Ij and we shall assume the 
natural condition that the volatility matrix a^{Ok,£k) is of full rank. We can then define 
the risk premium 

A'=(0fc,4) = a''{ek,£ky{aHOk,ikpHOk,iky)-'b''{0k,ik), 
which satisfies ()2.4p . 

Remark 2.3 One can write the dynamics of the assets model (12.2p -( i2?3]) - (12.5p as a jump- 
diffusion process under G. Let us define the G-predictable processes {bt)t>o and {crt)t>o 
valued respectively in M.'^ and W^^"^ by: 

n n 

^* = E^^-^i^'^'^'^'^)' ^* = ^'^n>^^^t{rk,Lk), (2.6) 

k=0 * k=0 * 

and the indexed G-predictable process 7, valued in M'', and defined by: 

n-l 

7t(£) = 5] If,. jH-Tk^LkJ). 
^ — ^ t- 

k=0 

Let us introduce the random measure fi{dt, di) associated to the jump times and marks 
{Tk,Lk), k = l,...,n, and given by: 

fii[0,t]xB) = ^l,,<tlL,GB, t>0,BeBiE). (2.7) 

k 

Then, the dynamics of the assets value process S is written under G as 

dSt = St*{btdt + atdWt+ [ -ftii)fi{dt,de)). (2.8) 

Je 

Notice that in the formulation ()2.8p . the process W is not in general a Brownian mo- 
tion under (P, G), but a semimartingale under the density hypothesis, which preserves the 
semimartingale property (also called (H') hypothesis in the progressive enlargement of fil- 
trations literature). We also mention that the random measure fi is not independent of W 
under the conditional density hypothesis. Thus, in general, we de not have a martingale 
representation theorem under (P, G) with respect to W and //. 
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In this market, a credit derivative of maturity T is modelled by a ^y-measurable random 
variable Ht, thus decomposed in the form: 

n 

Ht = Y.H^Tirk,Lk), (2.9) 

fc=0 

where H^{., .) is Tt ^ ^(A^) (g) S(£'^)-measurable, and represents the option payoff when 
k defaults occured before maturity T. 

The above model setup is quite general, and allows us to consider a large family of 
explicit examples. 

2.3 Examples 

Example 2.1 (Exogenous counterparty default) We consider a highly risky underly- 
ing name (e.g. Lehman Brothers) which may have an impact on many other names once the 
default occurs. One should take into consideration this counterparty risk for each asset in 
the investment portfolio, however, the risky name itself is not necessarily contained in the 
investment portfolio. A special case of this example containing one asset (without marks) 
has been considered in [5] , see also [131 H] ■ 

There is one default time r (n = 1), which may induce jumps in the price process S of 
the d-assets portfolio. The corresponding mark is given by a random vector L valued in E 
C [—1, oo)'^, representing the proportional jump size in the d-assets price. 

The assets price process is described by 

St = Silt<r + Sj:{T,L)lt>r, 

where 5" is the price process before default, governed by 

dS^ = S^ *{b'idt + a^dWt) 

and the indexed process S'^{e,l), (9,1) eR+xE, representing the price process after default 
at time 9 and with mark i, is given by 

dS}{e,e) = Sl{e,l)*{bj{9,l)dt + al{9,i)dWt), t>9, 
Sl{9,l) = S^e*{U + l). 

Here W is an m-dimensional (P, F)-Brownian motion, m > d, b^, are F-predictable 
bounded processes valued respectively in M.'^ and M*^^"*, and the indexed processes b^, 
lie in 7'if(M+,£'), and valued respectively in R'^ and R*^^™. 

Example 2.2 (Assets portfolio with multilateral counterparty risks) The defaults 
family and the assets family coincide, each underlying name subjected to the default risk of 
itself and to the counterparty default risks of the other names of the portfolio. The assets 
family is represented by a portfolio of defaultable bonds. Recall that a defaultable bond is 
a credit derivative which insures 1 euro to its buyer if no default occurs before the maturity, 
otherwise, the buyer of the bond receives a recovery rate at the default time. The recovery 
rate may be random, and so viewed in our model as a random mark at the default time. 

In this contagion risk model, the number of defaults times n is equal to the number d 
of defaultable bonds. We denote by P* the price process of the i-th defaultable bond of 
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maturity Tj, by Tj its default time, and Li its (random) recovery rate valued m E = [0, 1). 
The price process drops to Lj at the default time r^, and remains constant afterwards. 
Moreover, at the default times tj, j i (which are not necessarily ordered) of the other 
defaultable bonds, the price process P* has a jump, which may depend on tj and Lj. 
Actually, the jump size of P* will typically depend on Lj if the name i is the debt holder 
of name j. The assets portfolio price process S = (P^, . . . , P") has the decomposed form: 

n 

Pi = ^h,<t<r,^,P:'\T,,Lk,Lk), t>0, (2.10) 
k=0 

where Tk = (ti, . . . , fk) denotes the k first ordered times, = (ti, . . . , Lk) the correspond- 
ing index marks, i.e. fk = r^^, and = (L^^ , . . . , L^^.). The index F-adapted process 
P^'''{6k,tk,£k), for (0fc,tfc,£fc) € Afc X X E^, represents the price process of the i-th de- 
faultable bond, given that the k names {li, . . . ,6^) defaulted at times fk = Ok with the 
marks Lk = ik- Here, we denoted by = {{i-i, . . . , t^) € {1, . . . , n} : Lj ^ iji for j ^ j'}. 
When i € i.e. i = Lj for some j = 1,...,/e, then P'^'^{Ok, ik,t-k) = (-j, and 

otherwise it evolves according to the dynamics: 



dPi'\0k,Lk,ik) = Pi'\ek,Lk,ik)ibi'''iek,Lk,ik)dt + ai'\ek,Lk,ik)dWt), t>ek. 

Here W is an m-dimensional (P, F)-Brownian motion, m > n, and the indexed processes 
fji,k^ ^i,k Vwi'^ki^^ , E^), and valued respectively in and M^^"\ The jumps of the 

i-th defaultable bond are given by 



Pe^_^^^{(^k+i,i^k+i,£k+i) = Pg- {Ok,i^k,(-k){'^ + l0^^^{dk,i^k,f-k,i^k+iAk+i)), 



for 6*^+1 > Ok, and Lk+i G {1, . . . , ?i} \ {ii, . . . , tfc}, and we have T^^,^^ (^fc, tfc, 4, '-fc+i, 4+i) 

= -1 + 4+i/P*l'' (0fe,i.fe,4), meaning that Plf^UOk+i,Lk+iAk+i) = 4+1, when Lk+i = 

i. This model is compatible with several ones in the literature, see e.g. [2], [3], and we shall 
focus in the last section on this example for numerical illustrations in the case n = 2. 

Example 2.3 (Basket default swaps) A fe-th-to-default swap is a credit derivative con- 
tract, which provides to its buyer the protection against the A:-th default of the underlying 
name. The protection buyer pays a regular continuous premium p until the occurence of the 
/c-th default time or until the maturity T if there are less than k defaults before maturity. 
In return, the protection seller pays the loss 1 — Lk where Lk is the recovery rate if r^, is 
the A:-th default occuring before T, and zero otherwise. By considering that the available 
information consists in the ranked default times and the corresponding recovery rates, and 
assuming zero interest rate, the payoff of this contract can then be written in the form (12. 9p 
with: 

rrifa fi\ ( -p9k + {l-£k), iii>k, 

for Oi = {9i,...,9i) e Ai, = (4,...,^i) e E\ 
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3 The optimal investment problem 



3.1 Trading strategies and wealth process 

A trading strategy in the d-assets portfolio model described in section [2^2] is a G-predictable 
process vr, hence decomposed in the form 

n 

n = Vln^ vrf(rfc,Zfc), t > 0, (3.1) 
^ — ' *- 

fe=0 

where tt^ is an indexed process in 'P^{iS.f^, E^), and TT^{9k,'^k) is valued in A'' closed set of 
M.'^ containing the zero element, and representing the amount invested continuously in the 
d-assets in the /c-default scenario, given the past default events = 0^ and the marks at 
default Lk = £k, for {6k,£k) € x E^. Notice that in this modelling, we allow the space 
of strategies constraints to vary between default times. This means that the investor 
can update her portfolio constraint set based on the observation of the past default events, 
and this includes the typical case for defaultable bonds where the assets cannot be traded 
anymore after their own defaults. Notice that this framework is then more general than 
the standard formulation of stochastic control problem where the control set A is invariant 
in time. 

Remark 3.1 It is possible to formulate a more general framework for the modelling of 
portfolio constraints by considering that the set A'^ may depend on the past defaults and 
marks. More precisely, by introducing for any A; = 0, . . . , n, a closed set A^ c M'^ x x E'', 
s.t. (O,0fc,4) G A'' for ah (0^,4) G A^ x E^, and denoting by ^'=(0^,4) = {vr G M"^ : 
{7r,9k,£k) G a''}, the portfolio constraint is defined by the condition that the process 
7r'^(0fc,£fc) should be valued in A^{dk,£k). In the rest of this paper, and for simplicity of 
notations, we shall focus on the case where A'^ does not depend on the past defaults and 
marks, i.e. = x A^ x E^'^. 

In the sequel, we shall often identify the strategy tt with the family {Tr'^)k=o,...,n given 
in ()3.ip . and we require the integrability conditions: for all Of^ G A^, G E'^, 

\TT^i0k,ikybHOk,ik)\dt + / |7rf(6>fc,4)Vt'=(6'fc,4)|'dt < oo, a.s. (3.2) 

Jo 

where T < cxd is a fixed finite horizon time. Given a trading strategy vr = (vr'^)fc=o,...,n) the 
corresponding wealth process is defined by 

n 

Xt = Y.l^,X^{Tk,Lk), 0<t<T, (3.3) 

k=0 

where X^{8k,£k), ^fc S A^, 4 G E^, is an indexed process in Owi^k-, E^), representing 
the wealth controlled by ■K^{6k,(-k) in the price process S^{6k,(-k)^ given the past default 
events Tk = Ok and the marks at default Lk = £k- From the dynamics ()2.3p . and under 
(|3.2p . it is governed by 

dXf(0fc,4) = T:K^k,tk)'{h^t{^k,tk)dt + a\ek,tk)dWt), t>ek. (3.4) 
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Moreover, each default time induces a jump in the assets price process, and then also on 
the wealth process. From (j2.5p . it is given by 

Notice that the dynamics of the wealth process can be written as a jump-Ito controlled 
process under G by means of the random measure /i in p. 70 : 

dXt = 4{btdt + atdWt + [ jtii)tJ^{dt,di)). (3.5) 

Je 

3.2 Value functions and F-decomposition 

Let U be an exponential utility with risk aversion coefficient p > 0: 

U{x) = — exp(— px), x € M. 

We consider an investor with preferences described by the utility function U, who can trade 
in the d-assets portfolio following an admissible trading strategy vr € Aq to be defined below, 
associated to a wealth process X = X^''^ as in (|3.3|) with initial capital Xq- = x. Moreover, 
the investor has to deliver at maturity T an option of payoff Ht, a bounded ^T-measurable 
random variable, decomposed into the form (|2.9p . The optimal investment problem is then 
defined by 



V^{x) = sup E 



U{X^'^-Ht) . (3.6) 



Our main goal is to provide existence and characterization results of the value function 
, and of the optimal trading strategy vr (which does not depend on the initial wealth x 
from the exponential form of U) in the general assets framework described in the previous 
section. A first step is to define in a suitable way the set of admissible trading strategies. 

Definition 3.1 [Admissible trading strategies) 

For k = 0, . . . , n, A'^ denotes the set of indexed process tt^ in V^{/^k,E^), valued in 
satisfying (j3.'2|) . and such that 

• the family |c/(X^'(0fc, £fc)), r F — stopping time valued in [0fc,T]| is uniformly inte- 
grable, i.e U{X''{Ok,if^)) is of class (D), 



E 



Jl lEi-U)(x^iek,ik) + 7r'si0k,iky7^i0k,ikJ))m+ii£k,di)ds 



n — 1, 



< oo, when k < 



for all (0fc,£fe) e Afc(r) X E'', where we set Afc(T) = n [0,r]'=. We then denote by 
= {A^)k=o,...,n the set of admissible trading strategies vr = {7T^)k=o,...,n- 

As mentioned above, the indexed control sets A'' in which the trading strategies take 
values may vary after each default time. This nonstandard feature in control theory prevents 
a direct resolution to (|3.6p by dynamic programming or duality methods in the global 
filtration G relying on the dynamics (13. 5p of the controlled wealth process. Following the 
approach in [T3], we then provide a decomposition of the global optimization problem ()3.6p 
in terms of a family of optimization problems with respect to the default-free filtration 



10 



F. Under the density hypothesis (j2.ip . let us define a family of auxiliary processes € 
0^{/S.}^,E^), k = 0, ... ,n, which is related to the survival probability and is defined by 
recursive induction from a" = a, 

«t(0fc,4)= / / at+'(0fc,efc+i,4,4+i)d^fc+i%+i(4,f^4+i), (3.7) 

Jt JE 
for k = 0, . . . , n — 1, so that 



a- 



a 



t ' 



where dOf. = dOi . . . dOk, vid-^-k) = ^i('^^i) • • • %(4-i) c?4)- Given vr G A^, we denote 
by X^'^{6k,£k) the controlled process solution to ()3.4p and starting from x at 6*/;;. For 
simplicity of notations, we omitted the dependence of X^'^ in tt^. The value function to 
the global G-optimization problem (|3.6|) is then given in a backward induction from the 
F-optimization problems: 



= ess sup E 



ess sup E 

T 



C/(X^'"-//^)ar(0,^)|^e„J (3.8) 
f/(X^'^-//^)a^(0fc,4) + (3.9) 

A; + iK-j-iK-|-i 



for any x € M, A; = 0, . . . , n, {6k,ik) S ^k(T) x £^''. Here X*^'^ denotes wealth process in 
controlled by vr'^, and starting from x at time 6^. To alleviate notations, we omitted 



and often omit in the sequel in X'^'^, H^, tt^, 7'^ the dependence on (0fc,£fc), when there is 
no ambiguity. Notice that (dk^^k) appears in (|3.9p as a parameter index through X^'^, H^, 
it'', 7*^ and a^. On the other hand, 0^ appears also via 6^ as the initial time in ()3.9p . The 
interpretation of the relations ()3.8p - ()3.9p is the following. represents the value function 
of the optimal investment problem in the /c-default scenario, and the equality (|3.9p may 
be understood as a dynamic programming relation between two consecutive default times: 
on the /c-default scenario, with a wealth controlled process X^, either there are no other 
defaults before time T (which is measured by the survival density a'^), in which case, the 
investor receives the terminal gain U{X^-H!f), or there is a default at time r^+i, which 
occurs between 0^ and T, inducing a jump on X^ , and from which the maximal expected 
profit is V^^^ . Moreover, if there exists, for all k = 0, ...,n, some € attaining 
the essential supremum in ()3.8p - ()3.9p (independently of x), then the trading strategy vr = 
G Ag, is optimal for the initial investment problem (|3.6p . 



[J^'')k=o,. 



4 Backward recursive system of BSDEs 

In this section, we exploit the specific form of the exponential utility function U{x) in order 
to characterize by dynamic programming methods the solutions to the stochastic optimiza- 
tion problems ()3.8p - (|3.9p in terms of a recursive system of indexed Backward Stochastic 
Differential Equations (BSDEs) with respect to the filtration F assumed from now on to be 
generated by the m-dimensional Brownian motion W . 
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We use a verification approach in the fohowing sense. We first derive formahy the 
system of BSDEs associated to the F-stochastic control problems. The main step is then to 
obtain existence of a solution to these BSDEs, and prove that this BSDEs-solution provides 
indeed the solution to our optimal investment problem. 

Let us consider the starting problem (j3.8p of the backward induction. For fixed {0,£) 
€ A,i(T) X problem (13. 8|) is a classical exponential utility maximization in the market 
model S"'{6, £) starting from and with random endowment = H^+^ In ay- We recall 
briefly how to derive the corresponding BSDE. For t G [9n,T], € A^, let us introduce 
the set of controls coinciding with v until time t: 

An/ 4. ,,n\ ( „n r- A^- . „n ,,n ~l 

•A^{t,v ) = {vr € : vr.^t = i/.^tj, 
and define the dynamic version of (jS.Sp by considering the family of F- adapted processes: 



Vi^{x,e,£,v'^) = esssup U {X'^'"' - H^)\Tt , t > On, (4.1) 

so that Vf^ {x,6,£,i''^) = V^{x,0,£) for any i/"" € A^. From the dynamic programming 
principle, one should have the supermartingale property of {V[^{x,6,£,i'^),9n < t < T}, 
for any z/" G A^, and if an optimal control exists for (|4.ip . we should have the martingale 
property of {Vj"(x, 0, £, tt"). On < t < T} for some vr"" € A^- Moreover, from the exponential 
form of the utility function U and the additive form of the wealth process X" in p.4p , the 
value function process V"' should be in the form 

vr{x,e,£,u^) = u{xr-Yr{e,£)), 9n<t<T, 

for some indexed F-adapted process y" independent of z^", that we search in the form: dY/^ 
= —fpdt + Z'^dWt- Then, by using the above supermartingale and martingale property of 
the dynamic programming principle, and since Vf:{x,6,£,i'^) = U{x — H^) by (14. ip . we 
see that (y",Z") should satisfy the indexed BSDE: 

(En) yr(6>,£) = m(e,£) + -lnaTie,£) 

P 

+ / r{r,Z^,e,£)dr- [ Z^.dWr, en<t<T, 
Jt Jt 

and the generator /" is the indexed process in V^i^"^, A„, E"^) defined by 

r{t,z,e,£) = inf |g|^-ar(0,€)'7r|'-6"(0,£)V} (4.2) 

= -xi{e,£).z-^\\^{e,£)\^ + l ini \z + -\^{e,£)-cj^{e,£y^\\ 

where the second equality comes from (|2.4p . This quadratic BSDE is similar to the one 
considered in [15] or [6] in a default-free market. Next, consider the problems ()3.9p . and 
define similarly the dynamic version by considering the value function process: 

Vt''{x,ek,£k,u'') = esssup E [[/(X^'" - i/^(0fc, 4))4(^fe> 4) + (4.3) 



12 



for 9k < t < T, where A'^{t,u'') = {n'^ G : vr'x^ = iy%}, for u'^ G A'^, so that 
VQ^{x,0]^,tk,v^) = V^{x,6j^,lj^). The dynamic programming principle for (j4.3|) formally 
implies that the process 

for < t < T is a (P, F)-supermartingale for any G A^, and is a martingale for tt^ if 
it is an optimal control for (|4.3p . Again, from the exponential form of the utility function 
U , the additive form of the wealth process in (j3.4p . and by induction, we see that the 
value function process V'^ should be in the form 

v^\x,ei,,tk,y'') = c/(xf'"-y,'^(0fc,4)), ek<t<T, 

for some indexed F-adapted process Y'', independent of u^, that we search in the form: 
dY^^ = —ffdt + Z^dWt- By using the supermartingale and martingale properties of the 
dynamic programming principle for V'^, and since F^(x,0fc,4) = U{x - Hf^), with = 
a^, we see that {Y^, Z^) should satisfy the indexed BSDE: 



(Ek) Y^\ek,lk) = H^{ek,ik) + ^lna!^{Ok,£k) 

[■T pT 

+ f''{r,Y,\Z^,ek,ek)dr- Z^.dWr, 9k<t<T, 

with a generator /'^ defined by 
f''{t,y,z,ek,ik) = inf |||z-af(0fc,4)'^|'-fe'(0fc,4)'vr (4.4) 

+ ^[/(y)^f/(7r.7t'=(0fc,4,^) ->^/^+'(0fc,t,4,^))%+i(4,d^)} 



-Ai'(0fe,4).^-7^|Ai^X^fc>4)P 



2p' 

+ inf II z + -A^(0fe,4)-fTf(0fc,4)'vr' 

+ Uj{y) j^U{'K.^^{ek,lk,i)-Yt^+\ek,t,lkA))r,k+i{lk,di)], 



where the second equality comes from ([2 

The equations (Ek), A; = 0, . . . , n, define thus a recursive system of families of BSDEs, 
indexed by {0,£) G A„(T) x E'", and the rest of this section is devoted first to the well- 
posedness and existence of a solution to this system, and then to its uniqueness via a 
verification theorem relating the solution to the value functions (j4.ip . (j4.3p . 



4.1 Existence to the recursive system of indexed BSDEs 

The generators of our system of BSDEs do not satisfy the usual Lipschitz or quadratic 
growth assumptions. In particular, in addition to the growth condition in z for defined 
in ()4.4p . there is an exponential term in y via the utility function U{y), which prevents a 
direct application of known existence results in the literature for BSDEs. 
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Let us introduce some notations for sets of processes. We denote by S^[t,T] the set of 
F-adapted continuous processes Y which are essentiahy bounded on [t,T], i.e. ||y|| , 



esssup |^s(<^)| < oo, and by L?r,[t,T] the set of F-predictable processes Z s.t. 



{s,uj)£[t,T]xQ 

E[J^ \Zs\'^ds] < oo. For any k = 0, . . . ,n, we denote by 5^(Afc, E'') the set of indexed F- 
adapted continuous processes in Of{Ak,E^), which are essentially bounded, uniformly 
in their indices: 

rk\\ IW^kl 



We also denote by L^(Afc, E^) the set of indexed F-predictable processes in Vw{Ak, E^) 
such that 



E 



r \Z^{ek,ik)\''dt\ < oo, y{6kJk) eAk{T)x eK 



We make the following boundedness assumptions: 
(HB) 

(i) The risk premium is bounded uniformly w.r.t. its indices: there exists a constant C 
> such that for any A; = 0, . . . , n, (0fe,4) € Afc(T) x E^ , t e [6'fc,T], 

|Aj'(0fc,4)| < C, a.s. 

(ii) The indexed J-r-measurable random variables Hj, and Inoj. are bounded uniformly 
in their indices: there exists a constant C > such that for any /c = 0, . . . , n, (0^, 
€ Afc(r) X E\ 

|i^r(^fe,4)| + |ln4(^fc,4)| < C, a.s. 
We then state the existence result for the recursive system of BSDEs. 



Theorem 4.1 Under (HB), there exists a solution (y'', ^*'')fc=o,...,n € Jj5^(Afc,£' 

A:=0 

JjiY{Ak, E^) to the recursive system of indexed BSDEs (Ek), k = 0, . . . ,n. 

Proof. We prove the result by a backward induction on k = 0, . . . , n, and consider the 
property 

(Vk) there exists a solution G S^{Ak, E^) to (Ek). 

• For k = n. From the expression (j4.2p of the generatof there exists some positive 
constant C s.t. 

\r{t,z,6,£)\ < C7(|zp + |AJ'(6>,£)|2), V(t, z, 6>, £) G [0, T] x M'" x A„(r) x . 

Hence, under (HB), we can apply Theorem 2.3 in |12j for any fixed {0,i) G A„(T) x E^, 
and get the existence of a solution (y"(0, £)) G S^[en,T] x L2^[6'„,r]. Moreover, 
from Proposition 2.1 in [12], we have the following estimate 

|Ft"(6>,^)| < esssup(\HTie,i)\ + -\lnaTiO,i)\)+c[ \X^{e,i)\'^ds, en<t<T. 
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Under (HB), this implies that sup^^ ^^g^^^j-^^xE" ll^"(^'-^)lls°o[e„ t] ^ Fi^ia-lly, the mea- 
surabihty of y" and Z" with respect to {6,i) follows from the measurability of the coef- 
ficients H^, a^, and /" w.r.t. {6,i) (see Appendix C in [H]). The property (Vn) is then 
proved. 

• Fix k G {0, ... ,n — 1}, and suppose that (Pfc+i) is true, and denote by {Y^^^,Z^^^) € 
5^(Afc+i,-E''+^) xLp^,(Afc+i,£^''+^) a solution to (Ek+i). Since the indexed F-adapted pro- 
cess y'^+i is continuous, it is actually F-predictable, and so Y^^^ € Vvi^k+i^E'^'^^)- This 
implies that the map (t, w, 6>fc, 4+i) ^ Yl'^^{uj,ek,t,lk+i) is P(F) ® B{l^k) ® B{E^+^)- 
measurable. The generator is thus well defined in (14. 4p as an indexed process in 
P][r(M, R™, Afc, E'*^), and we shall prove that {Vk) holds true by proceeding in four steps, 
in order to overcome the technical difficulties coming from the exponential term in U{y) 
together with the quadratic condition in z for f^. 

Step 1: Approximating sequence. We truncate the term U{y) = — e~^^ when y goes to 
— oo, as well as the infimum, by considering the truncated generator 

f^it,y,z,ek,ik) = -x'iiek,ik)-z-^\xUek,ek)\^ + 

zp 

inf 11 z + -XUek,ik)-4i0k,£ky7r\ 

7r6A'=,|{crJ=)'7r|<Af ^ P 

-C/(max(-Ar,y)) / [7(^.7^ (0^, 4, ^) - >^t'+'(0fc, t, 4, (4, 

P J E ^ 

and introduce the corresponding family of approximated BSDEs with terminal data 
and generator f^: 



N- 

Y,''^^{Ok,£k) = F^(0fc,4) + -lna^(0fc,4) (4.5) 

P 



+ / f^{r,Y,^Z^''',0k,ek)dr- [ Z^^.dW,, ek<t<T. 
Jt Jt 

Under (HB)(i), there exists a constant C such that for all {Ok,£k) £ Afe(r) x E'^, 



fUt,y,z,GkAk) > -xH0k,ik).z-^\x1iek,ik)\^ 

2p 

> -C{l + \z\), (4.6) 

for all {t, y, z) G [Ok,T] xRxM™. Moreover, since € , and the process y^+i is essentially 
bounded, there exists some positive constant Cn (depending on N) s.t. for all {Ok,ik) G 
Afc(r) X E\ 



pi. 2 



+ Cn 



f^{t,y,z,ek,ek) < -x1{ek,ik).z-—\x'^{ek,ek)\^ + '^z + -x1{ek,£k) 

< Cn{1 + \z\^), (4.7) 

under (HB)(i). Hence, for any given {6k,£k) G Afe(T) x i?'^, we can apply Theorem 2.3 
in [ig, and obtain the existence of a solution {Y^'^ (Ok^k), Z^'^ {OkJk)) G S^[dk,T] x 
L'^[9k,T] to gSD. The measurability of (yfe.A',^^'^) w.r.t. its arguments (0fc,4) follows 
from the measurability of H^, a^, w.r.t. {6k,£k)- In the next steps, we prove the 
convergence of the sequence iY^'^ , Z^'^)^^ to a solution of (Ek). 
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Step 2: Lower hound for the approximating sequence. Define the generator function / by 

f{t,z,ek,£k) = -\Kek,e-k).z-^\\^t{Ok,t-k)?. 

2p 

Under (HB)(i), and for fixed (0^,4) G Afe(r) x E'', the function 0^,4) satisfies the 
usual Lipschitz (and a fortiori quadratic growth) condition in z, which imphes from Theorem 
2.3 in [II] that there exists {Y_^{9k,£k)^ Z^{Ok^(-k)) G S'^[ek,T] x lJ^[ek,T] sohition to 
the BSDE with terminal data H!^{6k,ik) + ^ lna|'(0fc,4), and generator /^(., ., 0^, 4). 
The solution {Y_^,Z_^) is measurable w.r.t. the arguments {Ok,ik), and from the uniform 
boundedness condition in (HB), and Proposition 2.1 in [12], we deduce that {Y^,Z^) G 
S'^{/\k,E^) X L^(Afc, £;'=). Moreover, we easily see under (HB)(i) that for any (0^,4) G 
Afc(T) X E^, f{.,ek,tk) satisfy Assumptions (H2) and (H3) of [E]. Since l%,ek,ik) < 
f^{-,Of^,£f^), we can apply comparison Theorem 2.6 in [12] to get the inequality: 

y/''^(0fc,4) > YUOkJk), 9k<t<T, a.s. (4.8) 

for all N, and (0fc,£fc) G Afc(r) x Since Y'' G 5~(Afc,S''), this implies that y'^'^ is 
uniformly lower bounded, and thus by (j4.5p . we see that for iV large enough, {Y^'^ , Z^'^) 
satisfies the indexed BSDE with terminal data H^, and with a generator where one can 
remove in the truncation in —N for U{y), i.e. 

>^/''^(0fc,4) = F^(0fc,4) + -ln4(^fc,4) (4.9) 

P 



+ f f^{r,Y,',Z^''',ekJk)dr- [ Z^'^'.dWr, 9k<t<T, 
Jt Jt 



with 



Zp 



+ inf |g z + -\HOk,ik) - (0fc,4)V ' 

+ -U{y) [ U{7r.^^iek,ek,i)-Y,''+\ek,t,ek,i))7],+M,d£)}, 
P Je ' 



Step 3: Monotonicity and uniform estimate of the approximating sequence. We cannot 

fk 

'AT' 



apply directly a comparison theorem for Y^'^ from the quadratic generators f%, and so 



we make an exponential change of variable by defining for any {Ok^^k) G Afc(T) x E^ , the 
pair of processes (y'='^(0fc, 4), ^"^'"^(^fc, 4)) G S^iOu^T] x L2,[0fc,r] by 

i;''^(0fe,4) = exp(py,'='^(0fc,4)) and zf'^(0,,4) = py^^(0fc,4)^f''^(0fc,4). 

A straightforward Ito's formula on ()4.9p shows that (y'^'^(0fc, 4); ■^^'^'^(^A:; 4)) is solution 
to the BSDE 

Yt^'^'iOk^tk) = a^T{^k,ik)exp{pH^{Ok,ek)) 

+ / /^(r,y,'='^,Z^^,0fc,4)fir - / Z^^^'.dWr, 9k<t<T, 
Jt Jt 
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where the generator is defined by 

f^it,y,z,0k,£k) = inf (0fc,4)'vr|' - p(A^X^fc,4)y + (^fc,4)' 

C/(7r.7f 4, ^) - 4, ^))r/,+, (4, . 

I J 



Fix {6k,ik) G ^fe X E''. Denote by g'^(Tr,t,y, z,6k,ik) the function inside the infimum 
defining i.e. /^(.) = inf^gAM{aJ^)'7r|<Ar S'atItt, ■)• Then, for all (t,y,y' , z, z' ,9k,£k) G 
[ek,T] xM? X (]R™)2 X Afc X we have 

\f^it,y,z,ek,ek)- fNit,y',z,ek,ek)\ 

< sup |5^(7r,t,y,z,0fc,4) -5'Ar(7r,t,y',2;',0fc,4)| 

7rGA'--,|{cr*^)'7r|<Ar 



< (^p'iV + piV|A^(0fe,4)|)|y-y'|+pA^k-^'l- 



Under (HB)(i), we then see that satisfies the standard Lipschitz condition in {y,z), 
uniformly in {t,oj). Since the sequence {f^)N is noninceasing, i.e. < f^, we obtain 

by standard comparison principle for BSDE that y'^'^+i < y^'^ , and so: 

y,'='^+^(0,.,4) < Y^^'^'iOk^k), ek<t<T, a.s. yN en, (4.io) 

for all {6]^,t]^) € Ay^ x . From the quadratic condition in z for /q in (|4.6|) - (|4.7|) . uniformly 
in {6k,£k)i and the a-priori estimate of Proposition 2.1 in jl2j . we deduce under (HB)(ii) 
that y'"'" G S^{Ak,E''). Together with (|4^ and (f4T0]l . this implies that there exists a 
positive constant M such that 

r'-'^lls-CA^i?^-) < M, ViVGN. (4.11) 

Step 4: Convergence of the approximating sequence. By using (j4.1ip in (|4.5p (or (j4.9p ). 
we see that (y'^.-'V^ satisfies the indexed BSDE with terminal data H^, and with 

generator given by 



fUt,y,z,0k,£k) = -x^iek,ik).z-^\x1iek,ik)\^ 

2p 



+ inf \l z + -X'l{9k,£k)-c7H0k,ikyTr 



:IK) 

+ -U{{-M)Vy) [ U{7r.j^{ek,£kJ)-Y,''^H0k,t,£kJ)H+^{£k,d£) 

P J E 

By the same arguments as for the generator there exists a constant Cm such that 

\f%{t,y,z,ek,£k)\ < CM{l + \zf), 

for all iV G N, {t,y,z) G [0,r] x M x M™, (0^,4) G x E''. Let us check that the 
nonincreasing sequence {f^)N converges uniformly on compact sets of {t,y,z) G [0, T] x 
M X M™ to f'' defined by 

f''{t,y,z,Ok,£k) = -XUOkJk)-z-^\XUOk,£k)\^ 

2p 

+ inf II z + -x^iek,£k)-a^iek,£ky7T ' 

+ -U{{-MVy) [ U{7r.j^{ek,£kJ)-Y,''+Hdk,tJkJ)H+d^k,d£)} 

V TP J 
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Indeed, notice that in the definition of f^, one may restrict in the infimum over vr in 
s.t. the function ^^(tt, .) inside the infimum bracket, i.e. 



g''{Tr,t,y,z,dk,ik) 



+ -f/((-MVy) / U{7r.j^{0k,£k,i)-Yt''+\ek,t,ek,i))v,+A^k,d£), 



is smaller than §^{11, .) for vr = 0. In other words, we have 



1 



,XU^k,ik)\'^ + ^ inf g''{TT,t,y,z,ek,ik) 



f\t,y,z,ek,ek) = 
where 

}C{t,y,z,9k,£k) = [tt eM.'^ : g''{7r,t,y,z,ek,£k) < 5*^(0, t, y, 2, 6>fc, 4)}- 

Since U is nonpositive, y'^+^ is essentially bounded, and under (HB)(i), there exists some 
positive constant C such that 



}Cit,y,z,ek,ik) C l^eM'^: z + iA^(0fc, 4) - ^'(^fc, 4)'^ < z + -x^iek,ek] 
I p p 



C {vTGM'^: |ai'=(0fe,4)'vr| <C(jz| + l)}, 



(4.12) 



for all {t,y, z,6k,£k) £ [0, T] x M x x x E''^. This shows that on any compact of 
{t,y,z) G [0,r] X M X M™, we have JC{t,y, z,ek,£k) C {vr : |(CT^^)'7r| < N} for N large 
enough, and so = f^, which obviously implies the convergence of {f^)N to /'^ locally 
uniformly on {t,y,z) £ [0, T] x M x R"*. We can then apply Proposition 2.4 in [12], which 
states that the sequence {Y^^^ {OkJk), Z'''^ {ek,£k))N converges in S^[9k,T] x L^[6'fc,T] 
to {Y^{dk,£k), Z''{Ok,£k)) solution to the BSDE with terminal data H^, and generator 
Z'^. The indexed processes {Y'^,Z'^) inherit from {y^^^ ^ Z'^'^) the measurability in the 
arguments {Ok,£k) S ^k ^ E^- Moreover, from ()4.1ip . we see that also satisfies the 
estimate 

ll^''ll5^°°{Afe,£;'=) < 

Hence, this implies that one can remove the truncation term —M in the BSDE with gener- 
ator satisfied by {Y'',Z''). Therefore, {Y^,Z^) G S^{Ak,E^) x l,'^{Ak,E'') is solution 
to (Ek), which ends the induction proof. □ 



4.2 BSDE characterization by verification theorem 

In this section, we show that a solution {Y^)k to the recursive system indexed BSDEs 
provides actually the solution to the optimal investment problem in terms of the value 
functions V^, /c = 0, . . . , n, in ()4.3p . As a byproduct, we get the uniqueness of this system 
of BSDEs, and a description of an optimal strategy by means of the solution to these 
BSDEs. 

Theorem 4.2 The value functions , k = 0, ...,n, defined in (j4.ip . (14. 3p . from the 
decomposition of the optimal investment problem (j3.6p . are given by 

V,'{x,0k,ek,i^'') = ?7(xf'"-yi'=(0fc,4)), Bk<t<T, (4.13) 



18 



for all X € M, (6>fe,4) € x E\ G A^, where {Y\ Z%=o,...,n e J{S^{I\k,E'') 



fc=0 



L^(Afc, is the solution to the recursive system of indexed BSDEs (Ek), k = 0, . . . ,n. 
Here, X^''^ denotes the wealth process in p.4p controlled by , and starting from x and 
Ok- Moreover, there exists an optimal trading strategy tt = (vr )k=o,...,n =0,...,n 
described by: 



(4.14) 



-u{Yt''{ek,£k)) f u{Tr.j^{ek,£kJ)-Y,''+\ek,tjkJ)H+A£k,d£)}, 



fork = 0,...,n-l, {Ok^ik) G Afc(r) x E^ t e [0fc,r], a.s., and 

7ri"(6>,£) G argmin ^^(61,^) + -Ar(6>,£) -^^"(e,^)^ 
/or k = n, {6,1) G A„(T) xE'',te [6'„,r], a.s. 

Proof. Step 1: We first prove that for all k = 0, . . . ,n, u'' e A^, [/(X^'^ - y^(0fc,4)) 

n 

> V\x,ek,lk,y^). Let (y^^'=)fc=o,...,„ G n'5c"(^'^'^^^) ^ L^(Afc,^^') be a solution to 

fc=0 

the system of BSDEs (Ek), = 0,... ,n. For any x G M, (0^,4) € Afc(T) x , G 
we apply Ito's formula to the process: 



£!l{x,Gk,lk,v^ 



+ 



f I u{x^'- + u^.^^{ek,ekJ)-Y,^+Hek,s,ek 



kJ))Vk+ii^k,di)ds, 



for k = 0,...,n, and ^"^{x, Ok, in, i^") := - Y/'XOn, in)) , foi k = n, and ek<t<T. 

From the dynamics of X^'^ and Y^, we immediately get 

d^^{x,ek,ik,iy'') 

= -u{x'i^^ -Y,\ek,ik))[{f^{t,Y^\zlek,ik)-gK^lt,Y^\zlek,ik))dt 

+ {a^,{ekAk)'y^ -z^).dWt 



where 
9t{'^^t,y,z,ek,ik 



+ -U{y) [ U{'K.^'^{ekAkA)-Yt^+\6k,t,lk,t))n,+,{h 

P JE 



for A; = 0, . . . ,71 - 1, and g^{TT,t,y, z, On, in) = f |2 - (0„, £„)'7r| - 6J'(0„, £.„)'7r for k 
= n. Since, by construction, f^{t,y,z,Ok,ik) = ^^iweA^ Qti'^^i'^^y^ ^^^k^ik), and recalling 
that U is nonpositive, this implies that the process {S,t{x,dk,ik,i^''), Ok < t < T}, is a 
local supermartingale. By considering a localizing F-stopping times sequence {pn)n valued 
in [dk , T] for , we have the inequality: 

^^':A,„{x,ek,ik,iy'')\J't] < ^^^p,S^,Ok,ik,i^''), ek<t<s<T. 
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Now, by definition 13.11 of tlie admissibility condition for , and since the processes , 
yfc+i g^pg essentially bounded, the sequence (CsAp„ ^fc' ^k, ^^))n is uniformly integrable for 
any s S [0^, T], and by the dominated convergence theorem, we obtain the supermartingale 
property of ^'^(x, Ok,£k, u^). Therefore, by writing the supermartingale property between t 
and T, and recalling that = + | Ina^, we obtain the inequalities 



Tt 



(4.15) 
(4.16) 



+ 



which hold true for any E A; = 0, . . . , ?i. 

Step 2: The process f^^ Zg{6k,£k)-dWs is a BMO -martingale, for any k = 0, . . . ,n, {6k, (-k) 
€ Ak(T)xE^. By applying Ito's formula to the process exjp{qYt^{Ok,£k)) with q > p between 
any stopping time r valued in [^fc,r] and T, and recalling the terminal data 1^ = = 
H!^ + i In a^, we get 



1 



q{q-p)¥. 



+ E 



exp [qY^{ek,ik))\Z'l{eMfdt 

exp (gy,'=(0fc,4))(/'(t,l^^^f,0fc,4)-||^fP)di 
exp {qH^{ek,£k)) - exp {qY^^{ek,£k)) 



Tr 



(4.17) 



By definition of in ()4.4p . and since 

yfc+i g S^{Ak+i,E''+^), there exists a constant C 
such that for all (t, y, z) € [0, T] x M x W^, 

fit,y,z,ek,£k) < ^\z\^-CU{y). 

Combining this last inequality with (j4.17p . we get 



1 



q{q-p)E 



exp {qYt^{ek,£k))\z!^{ek,£k)\^dt 



< qCE 



T 



e^v{{q-p)Yt\ek,£k))dt 



+ E 



Under (HB)(ii), and since Y^ G 5^(A,fc, this shows that there exists a constant C s.t. 



E 



< C, for any stopping time r valued in [0^, T], 



which is the required BMO-property. 

Step 3: Admissibility of tt''. Let us consider the functions g^, k = 0, . . . ,n, defined by 



g''{7r,t,uj,0k,£k) 



p 



-u{Yl'{ek,£k)) [ u{Tr.j^{ek,£k,i)-Y,''+H(^k,tJk,i))v,+Mk,de), 

P J E 
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for = 0,...,n- 1, and ^"(Tr, t, w, 0, = |Zf(0fc,4) + - cTt"(0, £)'7r| . Recall 

that the indexed F-adapted processes and Y^~^^ are continuous, hence F-predictable. 
Therefore, the map (tt, t, uj, 0k, 4) ^ g^{T^, t, u, 6k, 4) is B{W^) (g> V{¥) ® ^(A^) O B{E^)- 
measurable. Moreover, for k = 0, . . . ,n, {6k,(-k) ^ ^ , we recall from Remark 12.21 
that either a^{6k,ik) = and ^^{6k,(-k,^) = 0, in which case, the continuous function vr — )> 
gk{7r,t,uj,9k,£k) attains trivially its infimum for vr = 0, or a^{6k,0k) and '^^{Ok,ik,^) are 
in the form a^{ek,ik) = {a''{Ok,£k) 0), 7'=(0fc,4,^) = {j''{ek,ik,i) 0) for some fuU rank 
matrix a^{6k,(.k)- In this case, the infimum of g'^{'K, ■) over vr S is equal to the infimum 
over 7f G (ct'^)'^'' of function ^'^(Tf, .) where: 

/(7f,t,u;,0fc,4) = ^ Z,'=(0fc,4) + -Ai^0fc,4)-^ 

P Je 

for = 0,...,n- 1, and g''{TT,t,u,e,i) = \Zl'{ek,ik) + {6 , t) - . We clearly have 
gk{0,t,uj,ek,ik) < liminf^fe(7f,t,a;,0fc,4), 

|7r|— ^oo 

which shows that the continuous function vf — )■ g^{Tr,t,uj,6k,ik) attains its infimum over 
the closed set {a^YA^, and thus the function vr g^{T^, t, uo, Ok,(-k) attains its infimum over 
A^{6k,(-k)- By a classical measurable selection theorem (see e.g. [K]), one can then find 
for any k = {),... ,n, G Vw{^k, E^) s-t. 

%^{6k,(.k) e argmin g^{'K,t,6k,(-k), 6k<t<T,a.s. 

for ah (0fc,4) G Afe(r) X E''. Let us now check that the trading strategy tt — )fc=o,...,n 
is admissible in the sense of Definition 13.11 First, by writing that g^{Tt^,t,9k,ik) ^ 
g^{0,t,6k,ik), we get similarly as in (I4.12P the existence of some constant C s.t. 

(0fc,4)'vrt'=(0fc,4)| < C{1 + \Z^{ek,£k)\, Ok<t<T, a.s. (4.18) 
for all (0fc,4) G AkiT)xE'', k = 0,...,n. Since Z'' G L2^(Afc,E*^), and recalling (HB)(i), 



this shows that tt satisfies (13. 2p for all k = 0, . . . ,n. Let us denote by X the wealth 
process in ()3.4p controlled by vr*^, and starting from x at 6k- By definition of it^, we have 

/'=(^,y,^Z,^0fc,4) = g|^f-af(0fc,4)'vrtY-&'(^fc,4)'7rf (4.19) 



i[/(y,'=) / C/(7rf.7f(0fc,4,^)-n'+'(^fc>*,4,^))r?,+,(4,c?^)| 
P J E 



2' 

+ 

for A; = 0,...,n-l, and /"(t, Y^", , 6>, €) = - ^^"(e,^)'?^^"!^ - 6J'(6>,£)'7r4" for = n. 
Prom the forward dynamics of Y^, we can then write for all Ok < t < T: 

U(X^'^-Y,') = U{x-Y,\)£^(p{Z'-ia''yn'))Rl 

with 

S^PiZ' - ia^YT:'^)) = exp (p [\z^ - {a'J n',) .dWs - ^ T \Z^ - {a'Jn':\'ds) , 



21 



and 

^ Je^. Je 

for k = 0, ... ,n — 1, and Rf = 1. Now, from Step 2 and ()4.18p . the process fg^p{Z^ — 
((7*'')'7r'^).dVF is a i?MO— martingale and hence (see [IQ]), £^(^{Z^ — {a^)'Ti:^)^ is of class 
(D). Moreover, since U is nonpositive, we see that \R^\ < 1, and so \U{X^'^ — Y^)\ < 
U{x - YI^^).£^ {p{Z^ - {a^yTT^)^ , which shows that U{X^''' - Y^) is of class (D), and then 

also U{X^'^) since Y^ is essentially bounded. It remains to check that for all A; = 0, ... , n— 1, 



E 



< oo 



By definition of tt^' (which implies (j4.19p ). the process ^'^(rr, 0^, vr'^) defined in Step 1, is 
a local martingale. By considering a localizing F-stopping times sequence {pn)n valued in 
[9k, T] for this local martingale, we obtain: 

rTApn r 



E 



E 



U{X, 



k,x 
TAp„ 



yTAj-U{x-Y,\) 



< E 



Uix-Yl) 



since U is nonpositive. By Fatou's lemma, we get the required inequality, and this proves 
that 7c^ S Af, for any k = 0, . . . ,n, i.e. vr = {Tt^)k=o,...,n is admissible: tt G .Ac- 
Step 4: Since vr = {Tt^)k=o,....n is admissible, and recalling that the processes Y^ are 
essentially bounded, this implies that the local martingales ^^{x,Ok.,£k,Ti'^), k = 0, . . . ,n, 
are "true" martingales. Hence, the inequalities in ()4.15p - (|4.16p become equalities for u = 
TT, which yield: 



?7(xf'"-y,'=(0fc,4)) = E 

+ 



Uix'^''^ - HUOk,£k))a'T{ek,ik) 



(4.20) 
(4.21) 



C/(X^- + 7r^7sW - Y^'+\ek,s,£kJ)yh+Ah,di)ds 



for A; = 0, . . . , n, (Ok,£k) € Afc(r) x , t € [9k.,T], x € M. Let us prove the properties 
(|4.13p by backward induction on k = 0, . . . , n. For k = n, from the additive form of the 
wealth process X^'^ and the exponential form of the utility function U, we observe that for 
any t G [6'„,r], vr" G AY{t,i''^), the quantity 



E 



-UiX^'"") 



aT{0,i) 



Tt 



does not depend on the choice i^" G A^. By combining (|4.15p and (j4.20p . we then have 

-u{x:^''' -HJ^{e,£)) 



Jtie,£) := ess sup E 



< U[-Y,^{e,i)) = E 



-U{Xl 



Tt 



aT{9,£) Tt < Jr(0,^), 
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where we used in the last inequality the trivial fact that tt" G Af{t,7r^). This shows that 
U{ - Yf!^{0,£)) = Jp{e,e), and so ^^"(x, 0, z^") = UiX^'"" - ¥^{6,1)) for any z^" G A^, 
x € M, {9,i) € A„(T) X £^", which is the property ()4.13p at step k = n. Assume now that 
()4.13p holds true at step k + 1. Then, we observe similarly as above that for any t € [0^, T], 
tt'^ G Awit,v'^), the quantity 



E 



7-k,x\ 



E 



+ 



-U{X^'^) 

^k,x I k „,k 



^ f U{X^'^ + nt^^ii)-Y,'+\6k,s,eiJ)) 



-u{x; 



k,x\ 



is independent of the choice G Af. By combining (|4.16|) and (|4.21|) . we then have 



ess sup E 



+ 



T-k^xx 



-u{xl , 
< u{-Y^\ek,ik)) 

-U{X!^^'' - H^{0k,£k))aUOk,ik) 



Vk+iihidik+i)dOk+i\J't 



E 



-uix^n 

rk+l f yk,x I „,fc 



+ 



Ki: i^k:. + ^k,. (4+1), e.+i, 4+i] 



-U{Xl 



k,x\ 



Vk+iihidik+i)dOk+i\J='t 



< J'(0fc,4) 



where we used in the last inequality the trivial fact that tt^ G A^{t,Tr'^). This proves that 
U[ — Y^{6}^,ikj) = Ji{6f^,ik), and thus the property (|4.13p at step k. Notice that this 
representation of shows as a byproduct the uniqueness of the solution to the recursive 
system of BSDEs (Ek). Finally, the relations ()i:2T]) for t = 9k, together with (|i33]) . yield 



v''{x,e,i) 

V''{x,6k,£k) 



UiX'^''^ -H^)aT{e,£)\Te, 
U{x!^'''-H^)aUOk,ik) + 

yk+i (^M^ ^ ^k^^^ _^k^^^ {ik+i),Ok+i,ik+iK^, (4, d£k+i)dek+i\Te, 



T 

Ok JE 



which prove that tt = (vr'^)fc=o,...,n is an optimal trading strategy. □ 

Remark 4.1 We recall that in a default-free market Ito model for stock price S with risk 
premium A and volatility a, the optimal trading strategy (in amount) for an exponential 
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utility funtion U{x) = —e and option payoff Ht, is given by (see [15] or [6]) 



TTj € arg mm 

tt£A 



Zt + -Xt - iatYn 
P 



2 



where {Y, Z) is the solution to the BSDE: dYt = -f{t,Zt)dt + ZrdWt, Yt = Ht, f{t,z) 
= inf^gyi \z + ^At — ((Tt)'7rp. In our multiple defaults risk model, inducing jumps on the 

stock price, we see from (|4.14|) the influence of jumps in the optimal trading strategy 
within the /c-default scenario: there is a similar term involving the coefficients A*^ and 
corresponding to the default-free regime case, but the investor will take into account the 
possibility of a default and jump before the final horizon, and which is formalized by the 
additional term involving the jump size 7''. In particular, if 7'' is negative (in the one-asset 
case d = 1) meaning that there is a loss at default, then the infimum in ()4.14p will be 
achieved for a value tt^ smaller than the one without jumps. This means that when the 
investor knows that there will be a loss at default on the stock, he will invest less in this 
asset, which is intuitive. In the next section, we shall measure quantitatively this impact 
on a two-assets model with defaults. 



5 Applications and numerical illustrations 

For numerical illustrations, we consider a portfolio of two defaultable names, and denote by 
Ti and T2 their respective nonordered default times, assumed to be independent of F, so that 
their conditional density (w.r.t. F) is a deterministic function. We suppose that ti and T2 
are correlated via the Gumbel copula which is suitable to characterize heavy tail dependence 
and is often used for insurance portfolio. More precisely, we let P[ri > 9i,T2 > 92\J~t] = 
P[ri > 6*1, T2 > 62] = exp ( - {{aiOi)'^ + (026*2)'^)^/'') with 01,02 > and /3 > 1. In this 
model, each marginal default time Tj follows the exponential law with constant intensity 
Oj, i = 1,2, and the correlation between the two defaults is characterized by the constant 
parameter /3. The case (3 = 1 corresponds to the independence case and a larger value of /3 
implies a large linear correlation between the survival events /3'*(T) = corr(l|^j>2-}, 1|^2>T})- 
The default density of r = {ti,T2) is thus given by 

a^{e,,e2) = G{9,,e2)^!^0^n{9,, e2)^-'^(n(0i, 02) + /3 - 1) 

where G{6i,02) = P(ti > 0i,r2 > 62) = exp(— ti(0i, 02))- As explained in Section 2.1 and 
Remark 12. 1^ the case of ordered default times fi = min(Ti,T2), r2 = max(ri,r2) can be 
recovered by considering the marks (ti, L2) indicating the order of the defaults (n, r2). The 
density of (n, f2, /-i, /-2) is given by: 

a{0,i,j) = l{i=ij=2}"'^(6'i,6'2) + l{j=2j=i}a'^(6'2,6'i), 

for 6 = {61,62) € A2. Before any default, the price process = (5^'*^, 5^'^^) of the two 
names is governed by a two-dimensional Black-Scholes model with correlation: 

dS^ = Sf *{b^dt + a^dWt), 

where = (b^'^,b'^'^) is a constant vector in M^, is the constant matrix: 

_ ( ai.oyr^ ai'V \ 
^ " V a2'0 J' 



24 



with fj^''^ > 0, cr^'^ > 0, p G (—1,1), and W = {W^jW"^) is a two-dimensional Brownian 
motion. The associated risk premium is then given by A'^ = (A^''^, A^''^) with 



0" 



Once the name j defaults at time Tj, it drops to zero, but also incurs a constant relative 
jump (loss or gain) of size 7* G [— l,oo) on the other name i 7^ j. We denote by S^'^{9i) 
= S'^'^{9i,j) the price process of the survival name i after the first default due to name j 
^ i at time Tj = 61. We then have Sg^{di) = Sg^{l + 7*), and we assume that it follows a 
Black-Scholes model 

dsi'\ei) = si'\ei){b''^dt + a''^dBi), t>ei, 



with constants 6*'^ and o"*'^ > 0. Here is the Brownian motion: = \J\ — p^W^+pW"^, 
ij2 _ p^2_ Finally, after both defaults, the two names can not be traded anymore, i.e. S"^ 

= (Sl'2,S2'2) = 0. 

We consider the investment problem with utility function U{x) = — e~^^, without option 
payoff Ht = 0, without portfolio constraint, and solve the recursive system of BSDEs. Since 
all the coefficients of the assets price, and the density are deterministic, we notice that these 
BSDEs reduce actually to ordinary differential equations (ODEs). We start from the case 
n = 2 after the defaults of both names. The solution to the BSDE (En) for n = 2 is clearly 
degenerate: 

Y\e,iJ) = ilna(0,i,i), e = (61,92) e A2, i,j (£{1,2}, i^j. 

Let us denote by ¥'^''{61) = ¥'^{61, i), i = 1,2, the solution to the BSDE (El) after the 
first default due to name i. Notice that the auxiliary function q;^'*(0i) = a^(6i,i) defined 
in ()3.7p is given for i,j = 1, 2, i 7^ j, by 

aVie^) = r a{e,,92,i,j)de2 = ^^((a.^i)'' + 
Jt 9^ 

The function y^'* is then given by the solution to the ODE 

¥l^\9i) = -[/31na, + (/3 - l)ln0i + lln((a,0i)^ + {a,tf) - {{a,9^f + [a.tff'P] 

+ f^\s,¥}^\9^)ds, 

where 

1 ■ 1 fr^'i 2 (n IV'^ -12 1,, 

f^\t,y,9i) = ^ + inf i^l -a^'V + -e'P'^y-^^ a{9i,t,i, j)\ , 



for i,j G {1, 2}, i ^ j. For A; = 0, the survival probability is equal to 
aO = P[ti >r,r2 >r] = exp(-r(af +0^)^/^^), 
and the function ¥° to the BSDE (EO) is then given by the solution to the ODE 



¥,' = -^{4 + 4)^^+ f f{s,¥!)ds, (5.22) 
P Jt 
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where 



We perform numerical results to study notably the following parameters: the loss or 
gain at default, the default intensities, the correlation between the defaults and between 
the assets. We choose the parameters of assets as below and fix them to be the same in all 
our tests: 6^'° = fc^.o = o.02, o-^'° = cj^'O = 0.1, b^'^ = }?'^ = 0.01, a^'^ = a"^'^ = 0.2, p = l 
and T = 1. 

In Figure dl we present the optimal strategies vr = {tt^,t:'^) at the initial time before 
any default, for different values of loss or gain at default and of default intensity. In Figure 
dl we consider a symmetric case where the default intensities ai and 02, and the loss/gain 
7^ and 7^ are equal respectively, so are the same for vr^ and vr^. We choose the correlation 
parameter p = and (5 = 2. The optimal strategy is increasing with respect to 7, which 
means that one should invest less on the assets when there is a large loss of default. When 
7 = 1, the strategy converges to the Merton one, since in this case, the gain at default 
of the surviving name will recompense the total loss of the default one. Furthermore, the 
strategy is decreasing with respect to the default intensity. So when there is a higher risk of 
default, one should reduce her investment. In particular, if the default probability is high 
and the loss at default is large, then the investor should sell instead of buy the assets. Only 
when 7 becomes positive and the gain at default is large enough to recompense the default 
risks, she can choose to buy the asset again. 



Figure 1: Optimal strategy tt before any default vs Merton tt 
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1 



Figure [2] plots the evolution of the value function before default, i.e. t — )■ V^'{x) = 
—^-pi^-Yt') where is the solution of the equation (I5.22P and we have chosen x = 
in the test. We consider various values of 7 with the same parameters as above and let 
ai = 02 = 0.01, /? = 2. The survival correlation is equal to p^{T) = 0.5846. We observe a 
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larger value function when the gain at default (7 > 0) is larger. We also notice that the 
value function in a loss at default (7 < 0) situation outperforms the no- loss case (7 = 0), 
which means that one can take profit from a loss of the risky stock by a shortsale strategy. 



Figure 2: Value function 
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Figure [3] plots the evolution of the optimal investment strategy 7r(t) for t € [0,T], T = 
1, when there is a default event at time r = 0.6, the parameters being the same as in Figure 
[2] with two different levels of loss at default 7. We observe a jump of the trading strategy 
at the default time in both curves. When there is a larger loss at default, one should invest 
less from the beginning, however, after the default occurs, the trading strategies on the 
surviving firm become the same whatever the loss at default is. 



Table 1: Optimal strategies vr^ and tt"^ before any defaults with various 7 and default 
intensities. 



7 


-0.5 


-0.1 





0.5 


1 


Merton 


ai = 0.01, 02 = 0.1, /3 = 2 


= 0.2936 




7r2 


0.462 
-1.047 


1.659 
-0.709 


1.892 
-0.498 


2.621 
0.623 


2.832 
1.168 


2 
2 


fli = 0.1, 02 = 0.1, /3 = 2 


= 0.5736 






-0.353 
-0.353 


-0.210 
-0.210 


-0.147 
-0.147 


0.556 
0.556 


2 
2 


2 
2 


oi = 0.3, 02 = 0.1, 13 = 2 


p' = 0.4555 
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We present in Table [T] the optimal strategies at initial time before defaults for firms 
with different levels of default risks (ai ^02). We still suppose equal loss or gain at default 
(7^ = 7^). Similar as shown in Figured! when the default intensity ai of the first firm 
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Figure 3: Time evolution of the optimal strategy tt given a default 
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increases, one should reduce the investment on this firm. In the case of high default risks 
and loss at default, one should sell instead of buying the risky asset. However, the strategy 
on the second firm (the one with 02 = 0.1) will in general increase when its counterparty 
becomes more risky. 

Finally, we examine the impact of correlation parameters p and /9 on the trading strate- 
gies before any default. In the following test presented in Table [21 we fix ai = 0.01 and 
02 = 0.1. We observe that the correlation p between the assets will modify the benchmark 
Merton strategies. When p increases, the investment on the less risky asset goes in two 
directions: one should increase its quantity in the loss at default case and reduce it in the 
gain at default for the more risky asset, one should always reduce the investment. 

Concerning the parameter /3, when there is a larger /3 and hence a higher correlation be- 
tween the survival events, one should increase the investment in the less risky asset and 
decrease in the more risky one. 
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Table 2: Optimal strategy vr^ and tt^ with various p and (3. 
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p = 0.3, /3 = 2 
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